Orientational glass transition in Ceo- 
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We construct a model for qualitative description of the orientational glass transition in Ceo on 
the spin-glass theory basis. The physical origin of the frustration and the disorder is discussed. 
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As is established in a number of experiments, Cgo 
crystalizes in a face centered cubic (fee) structure. At 
ambient temperature the molecules rotate almost freely 
with centers on the fee lattice sites, so that the space 
group is FmZm. When the temperature decreases to 
T c w 260 K the first order orientational phase transition 
takes place: the sites of the initial fee lattice become di- 
vided between four simple cubic sublattices with its own 
preferable molecular orientation in each sublattice. The 
broken symmetry group is Pa3. 

Moreover the neutron-diffraction experiments [see, 
e.g., reviewsQQ have shown that the orientations in the 
ordered state are so that the electron-rich regions (the 
interpentagon double bonds) face the electron-deficient 
regions of the the neighboring Cgo molecule: the centers 
of pentagons or the centers of hexagons. 

Solid Ceo undergoes a well-known glass transition at 
T g = 90K when no orientational motion can usually 
be detected. 1 i 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 Each C 60 molecule has a sin- 
gle misorientation separated from the ground-state ori- 
entation by a large energy barrier E g . The energy dif- 
ference between these two orientations is U and E g » 
jj i i i 3 j 4 j io j i^ These two minima of the intermolecular an- 
gle dependent energy were shown to be much lower than 
the energies of other mutual orientations of the pair 
of molecules Below T g molecules orientations are 
forming glass and the occupation probability of misori- 
ented molecules is essentially frozen at value of ~ 17%. 
And this static orientational disorder persists down to 
very low temperatures. The thermal energy is too small 
compared with the energy threshold between the two 
states for further reorientation to be possible. 

In spite of the recent progress the present understand- 
ing of the intermolecular interactions in solid Ceo is still 
imperfect and no single model is able to describe correctly 
its whole physical properties>i2i±£ This paper discusses 
glassy orientational transition in Ceo- The glass state 
in Ceo usually was described purely phenomenologically 
as "freezing" of an ensemble of two-level systems. The 
step towards microscopic description of glassy orienta- 
tional transitions in ortho-para hydrogen was developed 
recently in Refs ll4lll5t the developed approach is use- 
ful because it takes carefully into account disorder and 
frustration. The symmetry of Ceo molecules essentially 
differs from the symmetry of hydrogen; the method used 
in Refs H3ll5l can not be directly used for glassy orienta- 
tional transitions in Cgo- We show how the approach of 



Refs ll4ll5l should be modified for Cgo case. 

Now we follow the consideration of the ReflTi! We 
consider the restricted number of allowed orientations 
instead of the continuous rotations. Let us take into ac- 
count in the energy only the orientations with pentagons, 
hexagons or double bonds directed towards 12 nearest 
neighbors in fee-lattice. The Ceo molecule is constructed 
in such a way that if 6 of its 12 pentagons (or 6 of its 20 
hexagons) face 6 nearest neighbor double bonds (P and H 
states) then 6 of its 30 interpentagon double bonds face 
the remaining nearest neighbors. Now the energy matrix 
elements can take only three values; Jq, the energy of 
the general mutual position, J p pentagon versus double 
bond, and Jh, hexagon versus double bond. Following 
RemT] and putting Jq = we obtain from Fig. (2b) of 
RefE] J P = -300K and J H = -110K. The large number 
of the nearest neighbor bonds allows for the mean-field 
description of the problem. 

We do not use the multipole expansion. So in the 
framework of our model calculations it is possible to build 
up the allowed functions using only the harmonics with 
1 = 6. We need only their transformation properties and 
effectively take into account higher order terms. 

Let us construct the functions Pj(w) and Hi(uj) ex- 
plicitly in terms of cubic harmonics K m = K§ iTn , m 

1.2 13 (see, e.g., Ref ll8|) . All functions Pi and Hi 

are the sums of K m , invariant under the icosahedral 
symmetry of the molecule (i.e., belonging to the Ai g 
representation of the icosahedral group Ih) if icosahe- 
drons are naturally oriented in one of 8 properly cho- 
sen coordinate systems. The states Pi(Hi) have 6 pen- 
tagons (hexagons) and 6 double bonds directed towards 
12 nearest neighbors along different [100] axes. Pi(o>) 
describes the molecule rotated from the standard orien- 
tation B (following Ref ll9|) about [111] axis through the 
angle 97.76125°. The angle for H x (w) is 37.76125°. The 
functions f^C^), ^(w) and Pi(aj) (or if 2, -H3, H4) are 
obtained from Pi(u>) (Hi) by subsequent counterclock- 
wise rotations of the molecule by 90° around z axes. 

If written in standard coordinate frame with cartesian 
axes along the cube sites these functions have the follow- 
ing explicit form: 

Pi(w) = apK^uj) + [3 p [K s (oj) + K 9 (u) + K w (uj)] 

+ 7p[A'ii(w) + A 12 H+A 13 M], (1) 
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P 2 (lu) = a P K x {u) + (3 P [-K & {u) + K 9 {lo) - K w (u)] 
+ 1 p[-K 11 {u) + K 12 {u J )-K 13 {lo)], (2) 

P 3 {u>) = a P K x {u) + (3 P [K 8 {u>) - K 9 {uj) - K 1q (cj)} 

+ jp[K u (u) - K l2 {u) - K l3 {u)}, (3) 

P 4 (w) = a P K x {u) + Pp[-K 8 (oj) - K 9 (u) + K w {lo)] 
+ -Yp[-K n (u) - K 12 (u) + K 13 (w)], (4) 

with a P = -0.38866; f3 P = 0.31486; j P = -0.42877. 
The functions -ffj(w) have the same form as Pi(oj) but 
with the coefficients a H = 0.46588; (3 H = 0.37740; ^ H = 
0.34432. The functions are normalized to unity. 

In the mean-held approximation by minimizing the ori- 
entational free energy one can obtain the nonlinear inte- 
gral equation for the one-particle orientational distribu- 
tion functions gi(w) for a molecule on ith sublattice of 
fee Ceo- In the neighborhood of the bifurcation point Tb 
we have a linearized system: 
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dio' [B(u},u}')h 2 (u}') + A(uj,uj')h 3 (u>') 
+D{uj,uj')hA{uj')) = Q, (5) 



where gi(w) = fjbhi{uS) + ...; fi = 1/T — and analo- 
gously for three other sublattices. 

Here A(ui,ui'), B(ui,ui'), D(uj,uj') are the sums of in- 
teractions over the nearest neighbors in the sublattices. 
For example, the sum in the plane perpendicular to the 
x axis can be written explicitly in the form 

D(u,u/) = 2{[(i\(w) +P 4 (w)) Jp + (FiH + H A (w))J H ] 
x [P 2 (u/) + P 3 (cj') + H 2 (uj') + H 3 {lu')] 
+ [P 2 (lo) + P 3 (uj) + H 2 {io) + H 3 (cj)] 

x [(P x (o/) + P 4 (u/))Jp + + H a (uj'))J h ] 

+ [(P 2 (lo) + P 3 (lu))J p + (H 2 (lo) + H 3 (lu))J h ] 

x [PiV) + P 4 (u/) + H x {u') + H 4 (u')]} 
+ + P 4 (oj) + H x {uj) + H 4 (u)} 

x [(P 2 (u/) + P 3 {J))J P + {H 2 (J) + H 3 {J))J H ]} . 

(6) 

We add the condition hi(u>) to transform one into an- 
other under the action of the cubic group rotation ele- 
ments which leave the fee lattice invariant. At the bi- 
furcation point Tb > nontrivial solution with broken 
symmetry appears, corresponding to the orientationally 
ordered phase. We have T b = 275K (T^ xp] m 260K) and 

hi (w) = aPi (w) + bHi (u) + cK x (w), (7) 

ft2(a>) = aP 3 M + 6fl3(w) + cKi(w), (8) 

h 3 {u) = aP A {u) + bHi(cu) + cKi{u), (9) 

hi(uj) = aP 2 {uj) + bH 2 {u) + cKi{u), (10) 



So we obtain^ the bifurcation temperature, the sym- 
metry of the ordered phase and the ratio p of the number 
of molecules with pentagon facing neighbor double bond 
near the phase transition in good agrement with the ex- 
perimental data: 



Pp 



PH 



a + b 
b 

a + b 



0.608, 
0.392. 



(12) 
(13) 



The experiments show [see, e.g., Refs|l||2j that the ratio 
of the number of molecules in those two states is about 
60 : 40 at the phase transition temperature, and increases 
when the temperature decreases and freezes at T g . 

It is convenient to reformulate our results. In the ori- 
entational ordered phase all sublattices are physically 
equivalent: they convert to each other when we rotate 
the crystal on 90° along the z-axis. The expressions for 
h 2l h 3 and hi formally coincide with h\(Cj) in their "own" 
system of coordinates that we label by tilde. Then we can 
write Eqs.JSJ) as a single equation as follows: 

hi(u>) + ^fj du'Efaw'Mu/) = 0, (14) 



where 



dw' [B(u,u')h 2 (uj') + A(w,w')h a (u')+ 

+D{u 1 u')h±(uj l )]= [ dti'E{&,ti')hi{Q'). (15) 



aa p + ban + c = 0. 



(11) 



The matrix elements of E are symmetrical. They de- 
pend on J p , Jh and the coefficients aprm, fip(H), 1p{h)- 
The E-matrix can be diagonalized: 

E{Q,Q') = ^^E r [K 8 (Q)+K 9 (uj)+K w {u>)]+ 

+E S [K n {u>) + K 12 (w) + K 13 (w)] | X 

E r [K 8 (u>') + K 9 (u>') + K 10 (Q')]+ 

+Es[K u (u>') + K 12 (u') + K 13 (Q')\ 

(16) 

The solution of the Eq. l|14|) is the following: 

h{Q) = r [K s (u) + K 9 (u) + K w (u)} + 

+ a[K 11 {u) + Ki 3 {u)+K 13 (u)]. (17) 

The coefficients in Eas. (ll6|) - (|17|l are determined at the 
bifurcation point Tb. 

Equations (fTI)- (|1 1|) and Eq. (|17|) give the analytical so- 
lution for the distribution functions near the bifurcation 
point. It is not difficult to find the solution of the basic 
nonlinear integral equations for the distribution functions 
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to our solution, 
* 0.7 from below 
> [Jh changes 
Two types of 



in wider range of temperatures. This solution preserves 
the shape of Eas. H7llll) and Ea. (CT7Jl with temperature de- 
pendent coefficients a(T), b(T) and c(T) [r(T) and s(T)}. 
But the relative ratio of the molecules number [that is de- 
termined by a(T)/b(T) and r{T)/s{T) ratio] weakly de- 
pend on the temperature that does not correspond to the 
experimental data. The two minima in the potential en- 
ergy [correspond to P and H orientations] are separated 
by rather high barrieriSiiiil, that is not taken into ac- 
count in our model, Eq.|0, so we introduce Jp = Jp(T) 
and Jh = Jh{T) to correct it4£ It is obvious that the 
pair interactions between molecules become less sensi- 
ble to their mutual orientations when the temperature 
increases because the libration increases. The theoreti- 
cal estimate^ shows that the libration amplitude about 
[111] direction increases nearly twice when the tempera- 
ture changes from T g to T c . 

We tried to fit experimental data-^ 
Ecis.fTjl- lfTTJl . It followed that when p P - 
then Tf, — > 150K from above and Jh - 
sign at this value of pp, but Jp < 
mutual molecular orientations are profitable above this 
temperature: a) the pentagon of one molecule versus the 
double bond of the neighbor molecule or b) the hexagon 
of one molecule versus the double bond of the neigh- 
bor molecule. Below this concentration, pp = 0.7, the 
frustration appears: (a)-interaction is profitable and (b)- 
interaction is not. Let us remind that in our model the 
presence of the barriers between the energy levels is taken 
into account indirectly: through the effective parameters 
Jp and Jh- 

Optical experiments show that the orientational dis- 
order at low temperatures leads to the nonhomogeneous 
lattice deformation, see RefsIH and refs. therein. The 
neutron scattering data show that at low enough tem- 
peratures a relatively large fraction of the intermolccular 
contacts were with double bonds pointing to pentagonal 
or hexagonal faces. Only a small fraction of molecules 
were oriented at random. But this picture never led to 
the low temperature phase, but only to partly frustrated 
structures^ It follows from the experimental data^ 
that there is small anomaly at T w 150K but it is still 
unclear wether this temperature corresponds to the ar- 
rest of the free uniaxial rotation or to a glass transition 
similar to that observed near 90K. Note that the discrete 
description of the orientational ordering, Eas.(|7fl- (|llfl . in- 
cludes as the answer the equiprobable distribution of the 
molecule orientations in the cubic lattice, Ki(iv), that 
corresponds to the rotations in cubic lattice [see, e.g., 
Refl20| in addition to H and P states. 

In summary at high temperatures the molecules travel 
slightly between different types of mutual orientations. 
When the temperature decreases the possibility to change 
the orientation for the molecule decreases because of large 
potential barrier between the lowest orientational states. 
Then the H-state of the molecule is rather rare phe- 
nomenon at low temperatures. The glassy state is ob- 
tained when the thermal energy is not sufficient to over- 



come the potential barrier that separates the two orien- 
tational configurations^. The ff -state of the molecule [6 
hexagons and 6 double bonds] is then profitable from the 
double bond side and is not from the hexagon side. The 
P-state of the molecule [6 pentagons and 6 double bonds] 
is then profitable from the double bond side and from the 
pentagon side. We see that the behavior of our system 
at low temperatures is similar to the behavior of the "di- 
lute" molecular crystals, e.g. ortho-para hydrogen>i2ii£i 
.ff-states here play the role of p — H 2 molecules in ortho- 
para mixture. 

The solution, Eq.(|T7I) [and Eas. (|7ITT|) ]. correctly de- 
scribes symmetrically oriented order at all temperatures. 
But the coefficients r and s (a, b and c) should be deter- 
mined from the experimental data for pp and pn- On 
this basis we suggest the following model for the glass 
description in Ceo- 

Let us consider a system of particles on lattice sites i, j 
with Hamiltonian 



(18) 



where Jij are quenched Gaussian interactions with zero 
mean, 



P(Jij) = 



1 



: exp 



2J 2 



with J = J/VN. 

U=c[Ki(Q) + Ki(Q) + Ki (u)] 
+d[K\ l {Co) + K{ 2 {u) + K{ 3 ((j)], 



(19) 



(20) 



where c and d depend on the ratio of the P- and in- 
states [d/c = —0.914 when pp — 0.83, i.e. at the pressure 
P = 0,T = T g ]. 

Using replica approach we can write the free energy 
averaged over disorder in the form: 



(F)j/NkT = lim - max/ 1 ]>> Q ) 2 +^ ]T {q aP ) 

a a>/3 



— lnTr^j/a} exp 



a>0 



(21) 



Here t = J/kT, Tr(. . .) = dip ft dcos{6){. . .). 

The saddle point conditions for the free energy give the 
glass and regular order parameters 



,a/3 



Tr 


U a U^exp 




Tr 


exp 







Tr 


U a exp (o 
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Tr 


exp 





(22) 
(23) 
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and the auxiliary order parameter 



P = 



Tr 


(F/ Q ) 2 exp (e 




Tr 


exp 









Here 



a a>/3 



(24) 



(25) 



In the replica symmetric (RS) approximation 23 the free 
energy l|21|l has the form: 



4 4 



G?Z ( Z ' \ , -i 

, exp In Tr 

V 2 



exp ( 6* 



Here 



(26) 



(27) 



The extremum conditions for the free energy l|26|) give 
the following equations for the glass and regular order 
parameters: 



q = dz 




(28) 
(29) 



and the auxiliary equation 



p = dz 



Here 



dz L 



Tr 


U 2 exp 


[>: 


)] 


Tr 


exp 







dz I z 
;ex P 7T 



2tt 



(30) 



(31) 



The replica symmetric solution is stable unless the repli- 
con mode energy A is nonzero. For our model we have: 

Ars = l~t 2 x 



dz L 



(32) 

The results of the calculations are illustrated in 
Figffl The order parameters do not go to zero because 



Tr 


U 2 exp fe 


)_ 




"Tr 


Fexp (§) 


Tr 


exp (e\ 




Tr 


exp (ej 




T/J 



FIG. 1: Order parameters and the heat capacity evolution 
with the temperature. Here d/c = —0.914, red solid curve is 
the heat capacity, dotted blue and dash green curves are the 
order parameters m and y/q respectively, violet dashed curve 
is the replicon mode A rcp i. The replica symmetry breaking 
occurs at the temperature To corresponding to A rcp i = 0; the 
glass transition temperature T g usually is very close to Tq. 



J duju 2k+1 (Lu) ^ 0, k = 0, 1, . . M The co-existence of the 
orientational ordered state (the order parameter m) and 
the glass (the order parameter q) is in agrement with the 
experimental datai 2 ^^ 

The replica symmetry breaking occurs at the temper- 
ature To corresponding to A rcp i = 0; the glass transition 
temperature, T gi usually is very close to To- The mildly 
sloping curve with the broad maximum for the heat ca- 
pacity qualitatively corresponds to the curve C° xp (T) ob- 
tained in Ref.3. 

It is well known that the temperature of orientational 
transition increases with pressure*^ When the pressure 
increases the number of P-states decreases, but the num- 
ber of iPstates rises. It is known that Jp = Jjj = J and 
pn = 0.5 when thepressure P w 0.25 Gpa and T ps 300K 
[see Fig. 13 in Ref 1]. Then we get from our model using 
the bifurcation condition: J = — 242K. If we assume that 
Jp and Jh depend linearly on P then we get: 

J P = -300(K) + 230(K/Gpa)P, (33) 
Jh = -110(K) - 530(K/Gpa)P. (34) 

Then if P = O.lGpa then we find T c = 281K and 
p P = 0.56; if P = 0.2Gpa then we find T c = 293K and 
p P = 0.52; if P = 0.3Gpa then we find T c = 309K and 
pp = 0.48. These results agree well with the experimen- 
tal datai for T c and p P . We believe that our approach 
describes well the orientational transition at small pres- 
sure. 
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When P -> 1.3Gpa it follows from Eqs.©-© that 
J P — > [it changes its sign here] and J# < 0. So our 
simple model for the orientational ordering becomes in- 
valid. Experiments showi*2^i& that at these pressures 
the number of P-states is very small. 

The orientational glass transition is hardly experimen- 
tally seen at P ~ 0.2 Gpa [pp w pg at all temperatures 1 ] 
and P > 1.5Gpa [pp <C 1] that does not contradict the 
above description of the glass transition. When p P fa p H 
the P- and iJ-states are both profitable and the analogy 



with diluted multipole systems becomes invalid. When 
the pressure is large enough pp becomes very small and 
there is no sense speaking about disorder and frustration. 
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